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1  BUBBLE  CLOUDS 


In  this  report,  we  survey  the  basic  ingredients  that  go  into  the  bubble 
cloud  hypothesis  for  the  enhanced  acoustic  backscatter  seen  at  high  enough 
frequency  and  wind  speed.  The  basic  picture  that  has  been  proposed  is 
that  spilling  waves  generate  foamy  water  which  is  then  subducted  down¬ 
ward,  modifying  the  local  sound  velocity.  One  proposed  mechanism  for  the 
downwelling  current  necessary  to  accomplish  the  required  subduction  is  the 
Langmuir  circulation  cell;  we  will  see  that  this  is  sufficient  but  other  mech¬ 
anisms  may  also  contribute.  Our  discussion  naturally  breaks  up  into  the 
hydrodynamics  issues  related  to  bubble  generation  and  subduction  and  to 
acoustic  issues  related  to  scattering  from  a  given  distribution  of  air  density 
below  the  surface. 

1.1  Spilling  Breakers 


It  is  well  known  that  as  the  amplitude  of  a  gravity  wave  increases,  the 
waveform  develops  a  relatively  sharp  peak.  If  the  peak  slope  gets  too  large, 
the  wave  will  break.  That  is,  the  water  will  be  unable  to  remain  in  a  laminar 
state  near  the  wave  crest  and  turbulence  will  set  in.  There  is  a  qualitative 
difference  between  a  plunging  breaker  in  which  a  double-valued  nature  of 
the  crest  is  quite  extreme  and  the  more  gentle  spilling  breaker  in  which  the 
instability  is  confined  to  near  the  crest.  This  difference  is  clearly  seen  in 
Figure  1-1.  There  is  a  general  expectation  that  spilling  breakers  are  by  far 
more  the  common  type  in  the  open  ocean. 
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Figure  1-la.  Spilling  breakers  in  the  N.  Atlantic,  in  wind  force  6  (from  Coles,  1967) 


Figure  1-lb.  A  deep  water  plunging  breaker  in  the  N.  Atlantic  (from  Coles,  1967). 


An  estimate  of  the  wave  steepness  at  which  breaking  will  occur  can  be 
obtained  by  considering  the  limiting  wave  amplitude/ wavelength  ratio.  To 
do  this,  we  can  use  a  boundary  integral  method  to  find  steady-state  gravity 
waves.  If  we  consider  two-dimensional  inviscid  flow,  we  can  represent  the 
velocity  as 

v  =  V  x  z9(x,  y) 

where  #  is  the  streamfunction,  x  represents  a  horizontal  direction  and  y  the 
vertical.  Since  V  x  v  =  0,  we  have 

V2tf  =  0 

the  only  source  for  ^  is  at  the  wave  surface,  due  to  the  fact  that  the  tan¬ 
gential  velocity  vt  (and  hence  the  normal  derivative  of  the  streamfunction) 
is  discontinuous  at  the  fluid-air  interface.  We  can  therefore  write 

V  =  j  Gvtds' 

where  G  is  the  Laplace  Green’s  function.  If  we  use  an  assumed  periodicity 
of  A  =  2x  to  define  a  length  scale,  the  proper  Green’s  function  is  [Kessler, 
Koplik  and  Levine,  1988] 

Go  =  — -~  log  (l  —  2cos(x  —  x')e~h~v'\  +  e~^v~y  I) 

47T  v  7 

-\{\y  -y'\-{y  -y')) 

for  infinite  depth  and 

G  =  G0  +  7-  log  (l  -  2cos(z  -  x/)e-|2W-v-v  l  +  e-2|2/f-y-v'1) 

4?r  v  7 

for  fluid  of  unperturbed  depth  H.  To  find  vt,  we  use  Bernoulli’s  law 

\  v2  +  gy  +  p!p  =  \B 
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where  B  is  Bernoulli’s  constant.  The  air  pressure  can  be  taken  to  equal 
zero  and  the  normal  velocity  equals  zero  in  the  moving  frame  of  the  wave. 
Hence  v?  —  B  —  2gy  (we  are  neglecting  surface  tension).  We  can  rescale  g 
to  1  by  proper  choice  of  time  scale.  The  equation  for  the  interface  follows 
from  substituting  vt  in  the  integral  equation  and  setting  =  0  everywhere 
on  the  air-fluid  interface.  This  equation  must  be  supplemented  by  the  area 
constraint  (i.e.  incompressibility). 

J  ya  dx'  =  0 

and  the  wave  phase  speed  c  is  determined  by 

~  —  cy  (y  — »  — oo)  (infinite  depth) 
ty(y  —  — H )  =  —cH  (finite  depth) . 

To  solve  these  equations,  we  parameterize  the  curve  by  a  finite  set  of 
points  equally  spaced  in  arclength  and  iterate  using  Newton’s  algorithm.  For 
small  amplitude  a  =  maxy  —  miny,  the  wave  is  purely  sinusoidal.  At  slightly 
larger  a,  we  recover  the  analytic  formula  found  by  Stokes 

c2  =  ^(l-M2*2). 

Finally,  at  a  ~  .44,  we  reach  the  wave  of  maximum  height  at  which  the 
solution  branch  ends  (see  Figure  1-2).  This  is  sometimes  referred  to  as  the 
120°  Stokes  solution.  Any  attempt  to  put  more  energy  into  this  wave  will 
invariably  cause  the  wave  to  break.  It  has  been  pointed  out  by  Banner  and 
Phillips  (1974)  that  surface  wind  drift  may  lower  the  steepness  at  which 
breaking  sets  in,  at  least  for  small  scale  waves. 

The  preceeding  discussion  has  been  for  a  fully  periodic  wavetrain.  In  the 
ocean,  a  more  typical  instance  is  the  existence  of  wave  groups  which  advance 
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Figure  1-2.  Wave  shape  immediately  before  breaking  (ka  =  .44). 
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at  the  group  velocity  Vg  =  |c.  Since  this  is  smaller  than  the  phase  speed, 
individual  waves  inside  the  group  advance  through  it,  grow  in  amplitude  and 
then  subside.  If  the  amplitude  at  the  group  center  becomes  of  order  of  the 
aforementioned  maximal  amplitude,  breaking  will  occur.  As  pointed  out  by 
Donelan,  Longuet-Higgins  and  Turner  (1972),  this  process  is  approximately 
periodic  in  time  with  period  equal  to  twice  the  wave  period.  This  offers  an 
explanation  for  some  anectodal  evidence  regarding  the  periodic  appearance 
of  whitecaps,  with  a  periodicity  that  depends  on  the  wave  speed.  If  verified, 
a  tendency  for  periodic  repeats  of  bubble  generation  with  a  rate  determined 
by  the  measurable  surface  waves  (whose  velocity  is  given  by  the  Doppler  shift 
of  the  Bragg  scattering)  could  be  a  useful  discrimination  method. 

Given  that  breaking  is  determined  by  having  the  waves  grow  to  a  max¬ 
imal  amplitude,  the  percentage  of  the  ocean  surface  covered  by  whitecaps 
should  depend  both  on  the  wind  velocity  U  (measured  at,  say,  10  meters 
height)  and  the  fetch  (Monahan  and  Monahan,  1985).  The  latter  depen¬ 
dence  has  not  always  been  looked  at  explicitly;  in  most  of  the  early  literature 
on  phenomenological  fits  to  whitecap  coverage,  the  assumed  dependence  was 
taken  to  be  (see  Monahan  and  O’Muircheartagh,  1986) 


w  =  BUa , 

with  U  in  m/sec.  This  led  to  estimates  of  a  ~  3.4  with  B  around  3.8  x  10-6. 
Little  systematic  dependence  on  atmospheric  stability  has  been  observed. 
Perhaps  the  most  sophisticated  treatment  of  the  phenomenology,  discussed 
in  Monahan  and  O’Muircheartagh,  predicts  the  comprehensive  formula  w  = 
w(U,  AT,  Tw,  d,  F)  where  A T  =  Twater  —  is  related  to  the  atmospheric 
stability,  Tw  is  the  actual  water  temperature,  F  is  the  fetch  and  d  is  the  wind 
duration. 
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The  basic  model  of  the  spilling  breaker  itself  is  due  to  Longuet-Higgins 
and  Turner  (1974).  A  cartoon  picture  of  the  flow  is  given  in  Figure  l-3a.  The 
basic  idea  is  that  the  turbulent  water  is  treated  as  a  distinct  fluid  which  slides 
down  the  forward  face  of  the  large  amplitude  wave.  This  fluid  is  less  dense 
than  the  underlying  “laminar”  fluid  since  the  turbulent  flow  incorporates  air 
bubbles  —  typical  estimates  of  the  density  of  such  self-aerated  flows  suggest 
that  the  density  on  a  30°  slope  might  be  from  70%  to  90%  that  of  pure 
water.  This  flow  is  complicated  by  the  shear  between  the  falling  turbulent 
flow  and  the  rising  (in  the  frame  moving  with  the  wave)  laminar  basement; 
this  shear  presumedly  tends  to  inject  water  into  the  turbulent  layer  and  also 
retards  the  motion  via  friction.  As  the  wave  amplitude  subsides,  the  gravity 
impact  forcing  driving  the  turbulent  flow  downward  decreases  and  the  foam 
is  carried  over  the  top  to  the  back  face.  The  slope  at  which  this  occurs  is 
in  some  ways  similar  to  the  angle  of  repose  for  granular  flows  at  which  an 
avalanche  will  be  suppressed  (as  we  go  to  9  <  On)  by  the  shear  friction  with 
the  underlying  solid  ground.  A  picture  of  the  entire  sequence,  taken  from 
Donelan  and  Pierson  (1987)  is  given  in  Figure  l-3b. 

It  is  fair  to  say  that  the  problem  of  predicting  the  flows,  air  densities, 
and  eventually  bubble  sizes  and  number  distribution  »n  a  spilling  breaker  is 
far  from  being  solved.  However,  it  is  probably  a  reasonable  guess  that  the 
typical  spilling  event  and  the  frequency  of  occurence  depends  most  strongly 
on  the  wind  and  on  the  fetch. 

1.2  Langmuir  Circulations 


Foamy  water  can  only  scatter  underwater  sound  effectively  if  certain 
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Figure  l-3a.  Sketch  showing  the  features  of  a  spilling  breaker  which  are  incorporated  in  the 
theoretical  model.  The  wave  is  moving  from  the  right  to  left  and  has  a  whitecap 
on  its  forward  face.  The  velocities  in  both  the  wave  and  whitecap  are  measured 
relative  to  the  wave  crest,  with  positive  direction  downwards. 
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Figure  l-3b.  Schematic  diagram,  with  the  vertical  scale  exaggerated,  through  the  centerline  of 
a  group  of  waves.  As  the  wave  on  the  ’eft  at  the  top  advances  from  t  =  f0  to  f  = 
f0  +  x,  it  steepens  and  forms  a  sharp  wk  -ge  (labeled  W)  at  t  =  t0  +  2x.  This  is 
followed  for  a  short  while  by  a  spilling  breaker  (b),  with  a  hydraulic  jump  0)  at  the 
toe  of  the  breaker,  as  at  t  =  ^  +  3x.  As  the  wave  decreases  in  height  on  progress¬ 
ing  through  the  group,  the  action  ceases  and  a  foam  patch  (F)  and  water  drops 
(D)  are  left  behind. 


conditions  are  met.  First,  there  must  be  significant  air  density  at  depths  re¬ 
moved  from  the  air/water  interface  which  acts  as  a  pressure  release  boundary. 
Secondly  there  must  be  some  non-trivial  horizontal  structure  to  the  bubble 
density  profile.  A  uniform  layer  of  aerated  water  will  just  present  a  lower 
“effective”  surface  and  have  no  appreciable  backscatter. 


A  typical  phenomenological  assumption  for  the  onset  of  acoustic  backscat¬ 
ter  is  fU2  >  104  where  /  is  the  frequency  in  Hz  and  U  the  windspeed  in 
m/sec.  Early  estimates  of  how  far  down  bubbles  would  be  expected  to  be 
observed  came  up  with  the  conclusion  that  the  bubble  hypothesis  could  not 
account  for  the  increased  backscatter;  roughly,  if  one  believes  that  a  signifi¬ 
cant  backscatter  can  occur  at  low  frequencies  (/  —  100Hz)  for  U  >  10  m/sec, 
this  requires  bubble  cloud  protrusions  of  order  5-10  meters  below  the  effective 
pressure  release  surface.  A  100  p  size  bubble,  assumed  to  suffer  low  Reynolds 
number  Stokes  drag,  will  have  a  rise  velocity  of  2.2  cm/sec:  clearly,  we  must 
have  a  vertical  downwelling  of  sufficient  magnitude  and  for  sufficient  dura¬ 
tion  that  within  a  bubble  lifetime,  bubbles  are  indeed  advected  10  meters 
downward.  For  bubbles  that  start  out  at,  say  100^,  an  estimate  for  the 
dissolution  time  is  (see  JASON  report  JSR-87-101) 

(isjk)  (d^h)  ■23seconds 

~1  minute  at  5  meters.  (Small  bubbles  have  shorter  dissolution  lifetimes 
but  a  smaller  rise  velocity  compensates  to  some  extent.)  For  the  required 
penetration,  we  thus  might  require 


(uj  —  2.2  cm/sec)  ~ 


5  —  10  meters 
(60)(2.5  minutes) 


which  means  a  downwelling  of  perhaps  10  cm/sec.  See  Thorpe  (1982)  for  a 
more  comprehensive  discussion  leading  to  a  similar  conclusion. 
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The  leading  candidate  for  providing  the  necessary  downwelling  is  the 
Langmuir  circulation  cell.  (There  may  also  be  occasional  downdrafts  due  to 
sudden  cooling  giving  rise  to  a  convective  instability,  but  these  seem  to  be 
rare.)  These  cells  are  associated  with  oft  observed  windrows,  long  parallel 
streaks  in  the  wind  direction  caused  by  convergence  zones  at  the  surface 
current.  It  has  been  noted  by  Thorpe  and  Hall  (1982)  that  waves  break  with 
equal  frequency  in  windrows  and  between  them  (i.e.  there  is  no  statistical 
correlation  between  Langmuir  circulation  patterns  and  wave  breaking)  and 
so  the  wave  breaking  can  be  thought  of  as  providing  a  uniformly  distributed 
source  for  the  Langmuir  current.  Assuming  that  we  have  some  empirical 
understanding  of  the  rate  of  foam  generation  (and  perhaps  an  idea  of  the 
bubble  size  distribution  if  there  is  no  one  specific  “typical”  size),  we  merely 
need  to  understand  the  causes  of,  and  patterns  in,  typical  circulation  cells. 
Unfortunately,  this  has  proven  quite  difficult. 

Before  continuing  our  discussion,  we  would  like  to  emphasize  that  to 
date  there  is  no  definitive  proof  that  Langmuir  cells  are  necessary.  It  is  con¬ 
ceivable  the  wave  breaking  by  itself  may  under  some  conditions  push  enough 
macrobubbles  downward  to  affect  acoustic  backscatter,  at  least  at  somewhat 
higher  frequencies.  Our  goal  is  to  outline  one  self-consistent  picture  of  what 
could  be  happening,  with  much  additional  effort  needed  to  confirm  or  inval¬ 
idate  this  scenario. 

A  review  of  the  possible  causes  of  Langmuir  circulations  was  presented 
by  Leibovich  (1983),  with  more  recent  measurements  by  Weller  et  al.  (1985) 
and  Smith  (1991).  A  schematic  picture  of  a  typical  flow  is  presented  in 
Figure  1-4;  the  typical  flow  velocities  are  5-10  cm/sec  (Figure  l-5a),  and 
typical  spacing  and  depths  believed  to  be  connected  to  the  mixed  layer  depth 
can  range  from  several  meters  to  several  hundred  meters.  For  the  typical 
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Figure  1-4.  Illustration  of  Langmuir  circulations  showing  surface  and  subsurface  motions 
(from  Leibovich,  1983). 
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Figure  l-5a.  Three  time  series  of  vertical  velocity  from  13  December. 
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sonar  of  interest  in  active  acoustics,  returns  from  bubble  clouds  in  individual 
descending  curtains  will  probably  not  be  resolved  and  one  will  instead  see  a 
composite  return. 

A  theoretical  model  which  adequately  describes  the  observations  of 
oceanic  Langmuir  cells  is  still  lacking.  The  most  likely  explanation  seems 
to  be  that  of  a  Stokes  drift  (caused  by  waves)  interacting  with  an  initially 
horizontally  uniform  current  (but  vertically  varying)  uc.  In  the  basic  hydro- 
dynamic  equations  of  motion,  there  is  a  term  of  the  form 

/  =  u.  x  (V  x  uc) 

.  duc 


if  both  the  Stokes  drift  and  the  current  shear  are  in  the  wind  direction  (say 
x).  This  is  like  a  gravitational  force,  pointing  downward  since  ^  <  0; 
it  of  course  can  be  statically  balanced  by  a  pressure  gradient.  If  however 
the  gradient  of  this  force  is  positive,  the  ocean  will  behave  like  an  unstably 
stratified  fluid  and  begin  to  convect  via  instability  growth. 

Now  ut  is  determined  by  the  square  of  the  wave  amplitude  a  times  the 
characteristic  wave  frequency  <7,  whereas  the  current  shear  will  be  determined 
by  wind  stress  to  be  uJ/1/7,  where  vj  is  some  effective  kinematic  viscosity 
and  u ,  is  a  friction  velocity.  Via  dimensional  analysis,  one  can  predict  a 
dimensionless  Langmuir  number 


(k  is  the  wavenumber)  which  governs  the  onset  of  the  instability  once  dissi¬ 
pation  is  taken  into  account.  Stratification  of  the  fluid  will,  of  course,  tend  to 
reduce  the  instability  and  gives  rise  to  the  notion  mentioned  above  that  the 
cells  only  extend  downward  within  the  mixed  layer.  Typical  windspeeds  for 
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the  onset  of  circulation  range  around  10  m/sec.  The  growth  of  a  Langmuir 
cell  in  an  ocean  experiment  after  the  wind  increased  from  8  to  13  m/sec  is 
shown  in  Figure  l-5b.  There  is  no  obvious  explanation  for  the  time  depen¬ 
dence  of  the  wave  vector  in  this  data. 


Most  of  the  measurements  of  Langmuir  circulations  have  been  accom¬ 
plished  by  high  frequency  sonar  scattering  from  below  the  surface  (actually 
scattering  from  the  bubbles!)  (see  Figure  1-5).  We  would  like  to  mention  the 
possibility  of  using  the  SAR  interferometry  technique  (Goldstein,  Barnett 
and  Zebker,  1989)  of  possibly  being  capable  of  resolving  the  surface  circu¬ 
lation  pattern  of  large  Langmuir  cells.  The  SAR  imaging  resolution  scale 
is  limited  mostly  by  velocity  bunching  by  the  ambient  surface  waves.  This 
“bunching”  is  due  to  the  mistaken  assignment  of  position  by  the  SAR  algo¬ 
rithm  due  to  the  motion  at  the  sea  surface.  Under  moderate  sea  states  one 
might  have  a  velocity  variance  ( v 2)  ~  (50cm/sec)2  which  translates,  at  say 
5  km  range,  100  m/sec,  to  a  resolution 


(.50) 


5  x  103 

100 


~  25  meters. 


The  interferometric  technique  consists  of  two  SAR  antennas  spaced  apart  by 
a  short  distance,  and  the  velocity  of  the  surface  current  is  found  by  comparing 
the  images  of  the  two  synthetic  apertures.  The  claimed  velocity  resolution  in 
a  recent  experiment  was  5-10  cm/sec,  roughly  the  same  order  as  the  circula¬ 
tion  currents  which  have  the  further  virtue  of  being  roughly  spatially  periodic 
and  therefore  easy  to  spot  by  looking  at  peaks  in  the  (spatial)  Fourier  trans¬ 
form. 
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Figure  l-5b.  Growth  of  Langmuir  cell  size.  Solid  line  is  40  m/hr  growth  rate  (from  Smith, 
submitted  1991):  wind  increases  from  8  m/sec  to  13  m/sec  at  7:20  A.M. 
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1.3  Acoustic  Scattering 


The  story  so  far  has  been  that  spilling  breakers  act  as  a  source  for 
foamy  water  which  occasionally  is  advected  downward  to  depths  of  5  to 
10  meters  by  Langmuir  circulation  cells.  As  far  as  the  acoustic  problem  is 
concerned,  what  we  really  need  is  a  plot  of  air  density  as  a  function  of  depth 
and  horizontal  position.  This  is  true  if  we  are  scattering  from  clouds  of 
bubbles  as  macroscopic  regions  with  changed  acoustic  velocities,  as  distinct 
from  any  resonant  contributions  associated  with  “macro-bubbles”.  Since 
this  seems  for  the  moment  to  be  the  most  likely  hypothesis,  we  will  limit  our 
discussion  to  this  case. 

Measurements  of  bubble  density  as  a  function  of  depth  is  again  accom¬ 
plished  by  high  frequency  sonar.  A  review  of  a  typical  experimental  setup 
and  typical  results  is  given  by  Thorpe  (1986).  In  Figure  l-6a  we  reproduce 
the  results  of  a  vertical  ranging  sonar  showing  typical  bubble  plumes  as  a 
function  of  time;  in  Figure  l-6b,  a  sidescan  sonar  records  the  horizontal  struc¬ 
ture  showing  the  connection  between  wavebreaking  events  and  bubble  plume 
formation.  In  this  study,  evidence  of  bubble  bands  were  found  above  wind 
speeds  of  7  m/sec,  roughly  consistent  with  the  expected  onset  of  Langmuir 
circulations.  Somewhat  surprisingly,  the  typical  spacing  between  bands  is 
only  of  order  5  m,  much  shorter  than  the  typical  large  cell  size.  This  might 
be  evidence  of  more  complex  circulation  patterns  (nested  cells,  e.g.)  or  of 
the  incompleteness  of  the  Langmuir  cell  explanation.  Soon  after  initial  for¬ 
mation,  the  bubbles  quickly  lost  any  velocity  imparted  by  the  wave  and  were 
merely  advected  by  oceanic  currents. 

In  a  similar  experiment,  Farmer  and  Vagle  (1989)  have  measured  the 
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Figure  l-6a.  Bubbles  observed  using  a  vertically  pointing  sonar.  The  sonograph  (top)  shows 
clouds  of  bubbles  below  the  surface.  Below  this  are  contours  of  log  Mvand  plots 
of  Mv  measured  at  six  levels  in  the  depth  range  bracketed  at  the  left  of  the 
sonograph.  The  wind  speed  is  shown  at  the  bottom.  The  wind  direction  was 
southwesterly,  the  fetch  exceeding  10  km,  and  the  air  temperature  was  1.75K 
below  the  water  temperature. 
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Figure  l-6b.  Sonograph  from  side-scan  sonar.  The  range  is  measured  along  the  surface  from 
a  position  immediately  above  the  sonar.  The  near-horizontal  streaks  are  due  to 
sound  reflected  from  bubble  douds.  The  wind  was  6.5ms  \  westerly.  Groups  of 
breaking  waves  can  be  seen  approaching  the  sonar  down  the  beam  in  the  295 
degree  direction. 


air  volume  fraction  as  a  function  of  depth;  their  graph  is  reproduced  here  in 
Figure  l-7a.  Although  the  volume  fraction  is  quite  small,  the  effect  on  sound 
speed,  given  via  the  index  of  refraction 

ni  =  i  23, 000^ 

l  +  z/10  meter 

<f>  =  air  volume  fraction 

is  still  capable  of  causing  scattering.  If,  following  Henyey  (1991),  we  model 
bubble  clouds  as  cylinders  with  typical  radii  of  1-2  meters,  this  depth  distribu¬ 
tion  can  be  converted  to  an  area  distribution  for  plumes  as  a  function  of  depth 
(Figure  l-7b).  The  prediction  of  large  plumes  every  1500  m2  means  that  for 
cells  of  size  100  m,  there  is  a  peak  downwelling  (or  perhaps  a  wave  breaking 
that  more  effectively  inserts  foamy  water  into  the  downwelling  flow)  every 
15  m  or  so  along  the  windrows.  More  careful  sonar  measurements  should  be 
able  to  selectively  search  for  large  plumes  (by  the  necessary  range  gating  in 
a  vertical  system)  to  see  if  this  is  at  all  reasonable. 

Given  bubble  plumes  determined  in  the  above  manner  from  the  mea¬ 
sured  volume  fraction  data,  Henyey  (1991)  has  given  a  convincing  demonstra¬ 
tion  that  the  enhanced  backscatter  (i.e.,  the  Chapman-Harris  (1962)  curve) 
could  be  accounted  for.  In  some  sense,  the  acoustic  calculation  is  by  far  the 
easiest  piece  of  the  puzzle;  with  the  exception  of  very  shallow  grazing  angle, 
multiple  scattering  effects  are  negligible  and  more  exotic  phenomena  (such 
as  localization  due  to  repeated  interactions  with  plumes)  highly  unlikely. 

In  some  more  recent  work,  Henyey  (private  communication)  has  pointed 
out  that  some  new  data  suggests  that  the  bubble  cylinder  radii  may  actually 
be  somewhat  larger  than  the  1-2  meters  originally  chosen.  Again,  various 
realistic  choices  of  larger  cylinders  still  give  fairly  consistent  answers.  One 
should  note  that  there  has  been  no  reported  evidence  of  strong  asymmetry  of 
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Figure  l-7a.  Air  volume  fraction  from  the  Fasinex  experiment,  extracted  from  the  results  of 
Farmer  and  Vagle  (1989).  This  data  constrains  the  microbubble  plume  model  at 
a  wind  speed  of  12  m/s. 


Figure  l-7b.  Model  prediction  for  the  plume  spedng.  An  experiment  which  can  resolve  5-m 
plumes  should  have  a  resolution  cell  no  longer  then  I0*m*. 
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should  note  that  there  has  been  no  reported  evidence  of  strong  asymmetry  of 
the  backscatter;  this  is  consistent  with  scattering  from  rare,  isolated  plumes 
but  would  possibly  contradict  a  scattering  mechanism  based  on  a  continu¬ 
ous  enhancement  of  air  volume  fraction  all  along  a  Langmuir  downwelling 
curtain.  At  the  present  level  of  sophistication,  all  one  can  really  say  is  that 
physically  reasonable  choices  of  clouds  of  microbubbles  consistent  with  sonar 
measurements  can  account  for  the  enhanced  backscatter. 

One  aspect  of  the  current  multi-step  approach  to  explaining  the  acoustic 
response  is  the  possible  sensitivity  of  the  result  to  an  almost  endless  set  of 
environmental  issues.  To  briefly  recap,  whitecap  coverage  will  depend  mostly 
on  wind  speed,  but  also  on  fetch  and  on  air  and  sea  temperatures.  Langmuir 
circulation  patterns  can  depend  on  swell  (which  causes  Stokes  currents)  and 
depth  of  the  mixed  layer,  in  addition  to  wind  and  wind  direction.  Any 
experimental  efforts  to  study  acoustic  scattering  must  be  cognizant  of  the 
need  to  carefully  determine  these  controlling  parameters. 
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